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THE EVOLUTION OF THE ACCRETION DISK AROUND 4U 1820-30 DURING A SUPERBURST 
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ABSTRACT 

Accretion from a disk onto a collapsed, relativistic star - a neutron star or black hole - is the mechanism widely 
believed to be responsible for the emission from compact X-ray binaries. Because of the extreme spatial resolution 
required, it is not yet possible to directly observe the evolution or dynamics of the inner parts of the accretion disk 
where general relativistic effects are dominant. Here, we use the bright X-ray emission from a superburst on 
the surface of the neutron star 4U 1820-30 as a spotlight to illuminate the disk surface. The X-rays cause iron 
atoms in the disk to fluoresce, allowing a determination of the ionization state, covering factor and inner radius 
of the disk over the course of the burst. The time-resolved spectral fitting shows that the inner region of the disk 
is disrupted by the burst, possibly being heated into a thicker, more tenuous flow, before recovering its previous 
form in ^ 1000 s. This marks the first instance that the evolution of the inner regions of an accretion disk has 
been observed in real-time. 

Subject headings: accretion, accretion disks — line: formation — 
neutron — X-rays: binaries — X-rays: bursts 



stars: individual (4U 1820-30) — stai's: 



1. INTRODUCTION 

Superbursts (e.g. Cornelisse et al. 2000; Strohmayer & 
Brown 2002; KuuUcers et al. 2002) are the more powerful and 
rarer relations of the Type I X-ray bursts, and are thought to 
be due to the nuclear burning of the ashes (principally carbon 
and rp-process nuclei) that remain after thermonuclear process- 
ing of the accreted light elements (Strohmayer & Brown 2002; 
Cumming 2003). One of the defining characteristics of a super- 
burst is its long decay time of a number of hours, three orders of 
magnitude greater than a typical Type I burst (Kuulkers 2003). 
This allows high signal-to-noise spectra to be accumulated at 
regular intervals throughout the course of the burst. 

4U 1820-30 is a compact low-mass X-ray binary with a pe- 
riod of only 11.4 minutes (Stella, Priedhorsky & White 1987), 
and lies within the globular cluster NGC 6624 at a distance of 
~ 7 kpc (Vacca, Lewin & van Paradijs 1986; Hesser & Shawl 
1985; Rich, Minniti & Liebert 1993). Type I X-ray bursts have 
been observed from this system for nearly 30 years (Grindlay 
& Gursky 1976). The accretion rate inferred from the persis- 
tent X-ray luminosity is 1.5 x 10^^ g (Cumming 2003). 
The first, and so far only, superburst seen from 4U 1820-30 was 
observed by the Rossi X-ray Timing Explorer (RXTE) on 1999 
September 9 (UT) during a scheduled monitoring observation. 
The initial analysis and interpretation of the data was presented 
by Strohmayer & Brown (2002). Time-resolved spectroscopy 
presented by these authors showed that an iron Ka line and 
absorption edge were present in the spectra, and that their prop- 
erties changed over time. A likely origin for these features is 
reprocessing of the X-rays by the accretion disk (Day & Done 
1991). This process, called X-ray reflection, is commonly used 
to explain similar features in the spectra of Seyfert galaxies 
(e.g. Pounds et al. 1990; George & Fabian 1991; Ballantyne, 
Vaughan & Fabian 2003). 

In this paper we present the results of fitting the time- 
resolved spectra of the superburst with detailed models of X-ray 
reflection. The next section describes the model calculations 
and details of the spectral fitting. The results are presented in 



Section 3 and discussed in Section 4. 

2. MODEL CALCULATIONS 

The X-ray reflection models were calculated with the code 
of Ross & Fabian (1993) (see also Ross, Fabian & Young 
1999; Ballantyne, Iwasawa & Fabian 2001; Ballantyne, Ross & 
Fabian 2001). A one-dimensional slab of uniform density gas 
is illuminated by a blackbody continuum (defined from 1 eV to 
100 keV) of temperature kT, and allowed to come into thermal 
and ionization balance. The emergent spectrum is calculated 
with a Fokker-Planck/diffusion method that naturally accounts 
for the energy redistribution of photons by Compton scatter- 
ing (Ross, Weaver & McCray 1978). The following ions are 
included in the calculations: C V - VII, N VI - VIII, O V 
- IX, Mg IX - XIII, Si XI - XV and Fe XVI - XXVII. Fluo- 
rescence and recombination lines from these ions are computed 
and transferred simultaneously with the continuum, resulting in 
a self-consistent prediction of the emitted spectrum. 

The compact orbit of this binary suggests that the donor is an 
evolved low-mass helium star (Rappaport et al. 1987), and thus 
the accretion disk material will be He-rich and H-poor (Cum- 
ming 2003). Therefore, we assumed the reflecting medium has 
the He, C, N, and O abundances of cool extreme helium stars 
from the tabulation of Pandey et al. (2001). The Mg, Si and Fe 
abundances were left at their solar values (Grevesse, Noels & 
Sauval 1996), but the effective metallicity is increased (by ^ 3) 
due to the negligible levels of hydrogen. The helium is assumed 
to be completely ionized everywhere in the reflecting region. 

The computed reflection spectra are parameterized by the 
ionization parameter ^ = AirFx/nne of the uradiated gas, 
where Fx is the flux of the incident X-rays and uhc = 
10^^ cm^'^ is the assumed number density of helium. Since the 
blackbody emission is also seen directly from the neutron star, 
the fitted spectrum is the sum of the blackbody and the reflected 
emission, weighted by a factor R, the reflection fraction. This 
parameter measures the relative importance of reflection in the 
spectrum, and can be related to the covering factor of the accre- 
tion disk. 
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Strohmayer & Brown (2002) describes the RXTE data and 
the process of accumulating the ^ 80 64 s spectra during the 
burst. The preburst spectrum was used as the background, 
and was thus subtracted from the data during spectral fit- 
ting (Strohmayer & Brown 2002). This persistent emission 
(amounting to a luminosity of 1.8 x 10"^^ erg s~^) is driven 
by the accretion onto the neutron star, so using this spectrum 
as a background impUcitly assumes that the burst will not ef- 
fect this underlying emission. As we see below, we cannot be 
sure that this does not happen; however, Strohmayer & Brown 
(2002) found httle difference in their spectral fits between the 
preburst background and the detector background. Since the 
flux of the superburst varies from 20 x to 3x greater than the 
preburst emission, then the shape of the background will not 
greatly effect the fit results. 

The spectral fitting was performed with xsPECv.ll.2.0bp 
(Amaud 1996). The model consisted of a grid of reflection 
spectra (read in using the at able command) and attenuation 
from neutral gas along the line of sight (using the wabs model), 
parameterized by the column density A^h- The Galactic absorp- 
tion toward 4U 1820-30 is 1.5 x lO^i cm'^ (Dickey & Lock- 
man 1990). The reflection spectra were blurred by the effects 
of relativistic motion within the Schwarzschild metric (Fabian 
et al. 1989). The inner radius of emission (rin) was allowed 
to vary, but the emissivity was fixed as r~^, where r is dis- 
tance along the accretion disk, and the outer radius was frozen 
at 200 GM/c^, where M is the mass of the neutron star. The 
inclination angle of the disk to the line of sight was fixed at 
30°. The fitted energy range was 3^0 keV, but was changed to 
3-15 keV 7000 s into the burst because of background dominat- 
ing at higher energies. The errorbars on the best fit parameters 
were computed using the 2a uncertainties for one parameter of 
interest. 

3. SPECTRAL FITTING RESULTS 

The results of fitting the reflection models to the superburst 
data are shown in Figure 1 . We find that an ionized reflector 
provides a good description to the spectra for most of the burst. 
During the first 1000 s the spectra are well fit with a large ^, 
indicating significant ionization of the accretion disk. Over the 
next ~ 500 s ^ decreases by more than an order of magnitude. 
For times between ~ 2500 and 3000 s the spectra cannot be fit 
by the reflection model, except for the last three points, where 
only very loose constraints could be placed on ^. Strohmayer & 
Brown (2002) showed the spectrum was hardening very rapidly 
during that period, which indicates that some rapid evolution of 
the system was ongoing, and is not well modeled by our equi- 
librium calculations. However, the recombination timescale for 
H-like Fe is irocomb ~ l/neaFc{T). The slowest timescale 
occurs at the highest temperature, which is 10^ K at the sur- 
face of the reflector, so that irecomb ^ 10^^ s, assuming the 
recombination coefficient of Verner & Ferland (1996). This 
is the same order of magnitude as the disk dynamical time at 
10 GM/c^. The recombination timescale of He II is ~ 0.1 s, 
much longer than the dynamical time, but faster than the 64 s 
spectra examined here. Therefore, the assumptions of ioniza- 
tion equiUbrium or fully stripped heUum cannot be the cause of 
the poor fits in this region of the lightcurve. A more fundamen- 
tal change in the X-ray source or reflecting geometry may be 
required. FinaUy, after the last Earth occultation passage, the 
system was revealed to have settled into a final decaying state, 
with ^ ~ 100 erg cm s~^, although one other non-equiUbrium 



event took place at ~ 6000 s. 

While the disk inclination angle was fixed at 30° for the spec- 
tral fits, there is a previous constraint of 35°-50° found by An- 
derson et al. (1997) from modulations in the ultraviolet flux 
(due to X-ray heating of the secondary). Relativistic disk fea- 
tures are also dependent on the accretion disk inclination angle 
(Fabian et al. 1989), so can provide a consistency check on the 
reflection hypothesis. In Figure 2 we plot the 68%, 90% and 
99% confidence contours for and the disk inclination angle 
5500 s into the superburst. This segment was chosen as it has 
a low reduced and strong reflection features. While the disk 
fit seems to prefer a lower inclination angle, it is consistent with 
the results of Anderson et al. (1997) at the 2a level. 

The observed evolution of ^ from highly to weakly ionized is 
to be expected since the observed flux decays over the course of 
the burst (Strohmayer & Brown 2002). More fundamental in- 
sights from this analysis arise from the variations of A^h, R and 
Tin over the first 2000 s of the superburst. The changes of both 
Tin, which measures the inner radius of the disk (as judged from 
the width of the iron line and edge), and R, the reflection frac- 
tion, strongly indicates that the accretion disk is significantly 
affected by the superburst out to a radius of ~ 100 GMjf?. In 
addition, the line-of-sight absorbing column also went through 
a similar variation, either as a result of the burst itself or the 
changes to the inner accretion disk. All three of these parame- 
ters follow a similar pattern where they increase (rin and A^h) or 
decrease (i?) before returning to close to their initial value. This 
behavior seems to be closely linked to the blackbody tempera- 
ture fcT, rather than the flux, which decreases monotonically 
during the burst. 

4. DISCUSSION 

A likely scenario to explain the observed behavior is that the 
accretion flow in the inner disk is altered by the superburst so 
that it loses the ability to produce reflection features. This could 
be accomplished if the scale-height of the disk, H, increased 
due to heating of the gas by the radiation field. This would de- 
crease the surface density of the disk, allowing it to be more 
easily and completely ionized by the burst radiation, removing 
any spectral features in the observable band, and effectively be- 
coming a mirror for the blackbody spectrum. As there would be 
no evidence for reprocessing from this innermost region, the es- 
timate of the inner radius would be from further out in the disk 
where iron can still produce an emission Une. The additional 
blackbody emission reflected from the ionized mirror, would 
dilute the reflection features, resulting in the low value of R. 
Alternatively, the lower surface density could reduce the opti- 
cal depth of the gas to the point where the inner disk becomes 
a very inefficient reflector, while only the outer disk would be 
optically thick enough for efficient reprocessing. 

The main evidence for this thermal effect is seen in the cor- 
relations between Tin, R and fcT. The Compton temperature of 
the gas is kTc « kT, and reached ~ 3 keV, about twice the 
maximum disk temperature (Shakura & Sunyaev 1973). Since 
H oc Csr"*/^ oc T^/^r''/^, where Cg is the local sound speed, 
then H. will increase by ~ 40% from its non-illuminated value. 
Changes to the disk surface density can take place over a vis- 
cous time which is ^ 1000 s at 100 GM jc^ if a 0.2, where 
a is the Shakura-Sunyaev viscosity parameter (Shakura & Sun- 
yaev 1973). This timescale compares favorably with the ob- 
served changes in R and rin. An increase in H may also ex- 
plain the evolution of A^h as the thicker disk may bring clumps 
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of disk material into the line-of-sight. Data which extend below 
the 3 keV limit of RXTE would be needed to test if the absorp- 
tion was caused by ionized or neutral material. Under this hot 
disk scenario, once the peak in kT is passed the disk material 
cooled to lower and lower temperatures, and thus contracted to 
its original size. Simultaneously, the amount of material avail- 
able to reflect increased, so R grew larger, but ^ dropped due to 
a combination of the increasing surface density and the decay- 
ing burst flux. 

The connection between the ionization parameter and 
the disk surface density S can be formulated as ^ « 
LxiTiHeH/r'^'E, where Lx is the luminosity of the burst, and 
mHo is the mass of a helium nucleus. Using the equations for 
a gas pressure dominated Shakura-Sunyaev disk with electron 
scattering opacity (Shakura & Sunyaev 1973), a peak burst lu- 
minosity of 3.4 X 10^^ erg s'^ (Strohmayer & Brown 2002), 
and a = 0.2, we find ^ w 100 erg cm s~^ for r = 15 GM/c^, 
which is consistent with the values of ^ found near the end of 
the burst (Fig. 1). From the observed decay in flux (Strohmayer 
& Brown 2002), the surface density of the disk must be ^ 100 
times less than the Shakura-Sunyaev value at the start of the 
superburst, but only ~ 6 times less near the end in order to ob- 
tain the necessary level of ionization. This indicates the surface 
density of the inner accretion disk was significantly lower at the 
peak of the burst but increased as the burst decayed, supporting 
the interpretation that the disk material was severely affected by 
the explosion. 

Another process that could lower the surface density at the 
inner edge of an accretion disk is the interaction between the 
radiation field and material flowing in the disk. Radiation can 
remove angular momentum from the disk material, causing the 
disk to empty inward, thus lowering the surface density. There 
exists only a few calculations of the interaction between an 
accretion disk and the radiation field from a luminous central 
source (Walker & Meszaros 1989; Walker 1992; MiUer & Lamb 
1993, 1996). These models predict that the radiation will ex- 
ert a significant torque on the accretion flow, which increases 
the accretion rate close to the central object. This effect may 
be enhanced in the case of 4U 1820-30 due to the power of 
the superburst (as compared to the typical Type I burst consid- 
ered in these models), as well as an interaction between the disk 
and the expanding shell that was ejected from the neutron star 
just prior to the superburst (Strohmayer & Brown 2002) - the 
radius-expansion phase (Lewin, Vacca & Basinska 1984). It is 



difficult to rigorously compare this idea to the observations be- 
cause detailed modeling of the disk-radiation interaction for a 
superburst or a radius-expansion burst has not been performed. 

By the end of the observation, the accretion disk has recov- 
ered from its disruption by the burst, and is reflective down to 
a radius of 10-20 GM/c^, which is curiously larger than the 
innermost stable circular orbit of 6 GM/c^. Since this is not 
an X-ray pulsar, the disk would not be truncated by the star's 
magnetosphere. Rather, this may be evidence for the disk- 
radiation interaction described above, and the material within 
10-20 GM/c^ is no longer dense enough to provide a reflec- 
tion signature. At these late times, the fits also require strong 
reflection, with i? > 1, indicating that the disk covers a large 
solid angle from the X-ray source, and thus is thin relative to 
the radius of the neutron star or has a non-negligible flare at 
large radii. The latter possibility is supported by the reflection 
from r > 100 GM/c^ observed early in the burst. If the inner 
disk expanded from the X-ray heating, then a significant disk 
flare might be necessary in order to allow the distant material 
to be irradiated by the blackbody emission. 

Dipping behavior is seen in the X-ray lightcurve of 4U 1820- 
30 between 6000-7500 s into the burst (see Fig. 7 in 
Strohmayer & Brown 2002). This kind of variability in low- 
mass X-ray binaries is often explained as variable absorption 
from a two-component ionized medium close to the neutron 
star (Frank, King & Lasota 1987). The dipping observed from 
4U 1820-30 could then be due to the reformation of the accre- 
tion disk, however, the timescale is incorrect, with the dipping 
seen ~ 2000 s after the probable recollapse of the disk. How- 
ever, there is an interesting correspondence between the times 
when dipping is observed and an increase in TVh derived from 
spectral fitting (Fig. 1). This may indicate that the dipping is 
caused by partially ionized material that is slowly recombining, 
but farther away from the accretion disk. An obvious source for 
such material is the gas that was cast off from the neutron star 
during the radius expansion phase. This material could be at a 
large distance from the system (10^^ cm if ejected at the escape 
velocity). 
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Fig . 1 . — Results of fitting the time-resolved superburst spectra of 4U 1 820-30 with ionized reflection models. The top panel shows the evolution of the reflection 
parameters: the ionization parameter 5, the reflection fraction, and inner radius of the disk. The lower panel describes the evolution of the absorption column, 
kT of the incident blackbody, and the reduced of the spectral fits. The reflection fraction, rin, and Nn seem to be closely related to the evolution of kT. Data 
were flt from 3^0 keV for the points to left of the dashed line, and from 3-15 keV for the points to the right. Time equals zero corresponds to about 100 s after 
the start of the burst. The hatched region denotes the interval where the spectra were hardening rapidly, and were not well described by the equilibrium reflection 
model. 
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Fig. 2. — Joint 68%, 90% and 99% confidence contours on rin and the disk inclination angle for a spectral fit 5500 s into the superburst. Modulations in the 
ultraviolet flux from 4U 1820-30 previously constrained the inclination angle to between 35° and 50° (Anderson et al. 1997), which is consistent with the results 
from disk reflection at the 2<t level. 



